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Abstract. 

Future cluster surveys will observe galaxy clusters numbering in the hundred 
thousands. We consider this work how these surveys can be used to constrain 
dark energy parameters: in particular, the equation of state parameter w and the 
non-adiabatic sound speed c 2 . We demonstrate that, in combination with Cosmic 
Microwave Background (CMB) observations from Planck, cluster surveys such as that 
in the ESA Euclid project will be able to determine a time-independent w with sub- 
percent precision. Likewise, if the dark energy sound horizon falls within the length 
scales probed by the cluster survey, then c 2 can be pinned down to within an order of 
magnitude. In the course of this work, we also investigate the process of dark energy 
virialisation in the presence of an arbitrary sound speed. We find that dark energy 
clustering and virialisation can lead to dark energy contributing to the total cluster 
mass at approximately the 0.1% level at maximum. 
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1. Introduction 

Current observational evidence for the accelerated expansion of the universe is 
compelling [1-3]. This accelerated expansion is generally attributed to the presence 
of a "dark energy" component in the universe's energy budget, a component with a 
strong negative pressure which dominates gravitational physics on large scales. The 
nature of this dark energy is, however, still a secret the universe keeps to itself — at least 
for now — and much effort is being invested to measure the properties of this dark energy 
in order that we may reveal its secret one day. 

Because of its unknown nature, there is no shortage of explanations for what 
this dark energy may be [4]. Therefore, in order to distinguish between the various 
explanations, we must appeal to a range of observational tests, each probing a different 
aspect of dark energy dynamics. The distribution of galaxy clusters in the universe, as 
quantified by the cluster mass function, is one such observational test [5-8]. 

The cluster mass function (or halo mass function) is sensitive to dark energy 
properties primarily through their effects on the growth function. Firstly, the dark 
energy equation of state w = Pq/pq, where Pq and pQ denote respectively the mean 
pressure and the mean energy density, affects the universal expansion rate. This in 
turn alters the overall growth rate of structures, and can be probed by way of the 
cluster abundance as a function of redshift. Secondly, if the dark energy is associated 
with a characteristic scale (such as a sound horizon), then the growth function will also 
inherit this scale-dependence. Indeed, as demonstrated in [8-11], a non-adiabatic dark 
energy sound speed c 2 s differing from unity leads to the clustering of dark energy on 
scales larger than the sound horizon. This in turn introduces a scale-dependence in the 
clustering of matter, which, as we shall show, ultimately translates into a distinctive 
mass-dependence in the cluster abundance. Note that in the case of a vanishing sound 
speed, dark energy is comoving with dark matter and clusters on all scales, a case that 
has been studied in e.g., [12, 13]. 

The purpose of this work is therefore two-fold. Firstly, we wish to estimate as best 
we can the dependence of the cluster mass function on the dark energy equation of 
state and non-adiabatic sound speed. To this end, we make use of the Press-Schechter 
formalism [14], which takes as an input the linear matter power spectrum and the 
linear threshold density from the spherical collapse model. Linear power spectra for 
nonstandard values of w and c 2 . are readily calculable using a Boltzmann code such 
as CLASS [15]. Likewise, a method for estimating the (cluster- mass dependent) linear 
threshold density in the case of c 2 ^ 0, 1 was previously established in [9]. What remains 
to be determined is the fractional contribution of clustering dark energy to the total halo 
mass, and for this purpose we must address the issue of virialisation in the presence of 
clustering dark energy. Virialisation of a combined matter and dark energy system has 
been studied in, e.g., [16], in the limit c 2 s = in the context of the spherical collapse 
model. In this work, we propose a simple method to track the virialisation process in a 
spherical collapse for an arbitrary dark energy sound speed. 



Confronting the sound speed of dark energy with future cluster surveys 



3 



Secondly, having computed the cluster abundance as a function of w and cj., we 
wish to investigate the sensitivities of future cluster surveys to these parameters. In 
particular, the Euclid project expects to observe some half a million clusters through 
weak gravitational lensing [17]. Using the appropriate survey specifications, we perform 
a parameter sensitivity forecast for Euclid based on the Fisher matrix approach. 

The paper is organised as follows. In section 2 we review the spherical collapse 
model and the issue of virialisation in dark energy cosmologies. In section 3 we present 
the cluster mass function as functions of the dark energy parameters w and c 2 s . The 
sensitivities of Euclid to these parameters are investigated in section 4. Section 5 
contains our conclusions. 



2. Virialisation in the spherical collapse model 

2.1. The spherical top hat and equations of motion for the matter component 

In the spherical collapse model, a spherically symmetric overdense region is assumed 
to sit on top of an otherwise uniform background matter density field. The overdense 
region has a physical radius Ri = -R(rj) at the initial (conformal) time Tj, and a uniform 
initial energy density 

Pmfr) =^0(1 + (2.1) 

where p m (r) denotes the energy density of the background matter field. This is our 
spherical "top-hat" perturbation, and the mass contained within the top-hat region is 
given by 

Air Att 
M halo = -jPm(n)(l + C,M = yMr )(l + 5%)Xl (2.2) 

where r denotes the present time, and we have defined X = R/a as the comoving 
radius of the top hat. 

The evolution of the physical top-hat radius R with respect to cosmic time t is 
described by the familiar equation of motion 

1 d 2 R AnG , , h , h th 



rW—— «+<«+ 3 ^ < 2 - 3 > 

where we have incorporated in the equation a dark energy component denoted by the 
subscript Q. Like the matter component, the dark energy component is taken to be 
uniform inside the top hat region defined by the radius R. Equation (2.3) can be 
equivalently expressed as an equation of motion for for the comoving top-hat radius X 
with respect to conformal time r, 

| + ?4 = -^r a2 W£ + PQ(1 + 3c2)#], (2.4) 
where • = d/dr, T-L = aH is the conformal Hubble parameter, and 
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is the square of dark energy sound speed defined in the rest frame of the dark energy 
fluid. Note that in identifying c s in equations (2.4) and (2.5) as the rest frame sound 
speed, we have implicitly assumed that we are dealing only with length scales much 
smaller than the Hubble length, and that the relative velocity between the dark energy 
and the dark matter fluids is much smaller than the speed of light. We also assume c 2 s 
to be constant in time and space. 

Since the total mass of matter inside the top hat M halo = (47r/3)p^-R 3 is conserved, 
the top-hat matter density contrast 5^ can be easily expressed as a function of the 
top-hat radius, 



rth/ \ _ Pm( T ) 
Pm{T) 

= a + d) 



a(n) R{t) 



- 1 = (i + o 



X(r) 



3 



- 1. (2.6) 



In the absence of dark energy perturbations, equations (2.4) and (2.6) form a closed 
system for the matter component. 

2.2. Equations of motion for the dark energy component 

On sub-horizon scales, the pseudo-Newtonian approach [18] applies, so that the (Fourier 
space) equations of motion for the dark energy density perturbations read [9] 

5 Q + m{c 2 s - w)5 Q + (1 + w)9 Q = 0, (2.7) 

1.2 2 

§ Q + (1 - 3c 2 )W Q - y^~ 5 Q + ^Ga 2 [pJ m + p Q d Q } = 0, 

where the top-hat matter density contrast is given in Fourier space by 

Att 

5 m (k,r) = Y [(l + 5^ t )X?-X*]W(kX), (2.8) 
with the filter function 

W(kX) = —— \sm(kX) - kX cos(kX)] . (2.9) 
[kX ) 6 

The dark energy density contrast ^{k, r) is then related to its top-hat average <5q(t), 
as appearing in equation (2.4), via 

5$(t) = ^j dk k 2 W(kX)~5^(k, t). (2.10) 

We refer the reader to reference [9] for a more detailed discussion of this formulation, but 
note that in equation (2.7) we have only kept terms up to linear order in the perturbed 
quantities 5q and 9q. This is a reasonable approximation for 5q not exceeding unity, a 
condition that is generally satisfied during most of the spherical collapse process except 
the final moment of collapse [9]. 

Equations (2.4) and (2.7) have been solved in [9] for a range of top-hat masses to 
determine the linear threshold density 5 C as a function of the halo mass. In particular, 
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figure 5 in the said paper shows a distinct step-like feature in S C (M), with the step 
occurring in the vicinity of the "Jeans mass" 
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Mj (a) = —p m (a)^-^) , (2.11) 

where Aj = 2n/kj and kj = T-L/c s are the Jeans' length and wavenumber respectively 
associated with dark energy clustering. 

2.3. Dark energy contribution to the cluster mass 

So far we have not considered the possibility that clustering dark energy might also 
contribute to the mass of the collapsed object. Following [12], we define the dark energy 
contribution to the total mass as 

M Q = fR(T) 3 5p%(r), (2.12) 

with 

H h to = PQ to S Q to = PQ to jf to r )' ( 2 - 13 ) 

where we have smoothed the dark energy density contrast over the spherical top 
hat. In defining Mq in this manner, we have implicitly assumed that only the dark 
energy perturbations contribute to the total mass of the bound object, and that once 
gravitationally bound, the dark energy perturbations behave like non-relativistic matter. 
It then follows that 

M tot EE M halo + M Q (2.14) 

gives the total mass of the system. 



2.4- Virialisation 

In the spherical collapse model, a gravitationally-bound halo is said to have formed in 
the top-hat region at the instant at which the virial theorem is first satisfied. After this 
time, the spherical region is taken to remain in a fixed configuration, with a constant 
radius i? v i r and mass M v j r (even though formally the equations of motion dictate that 
the region should eventually collapse to a point of infinite density). 

When only dark matter clusters and virialises, mass is conserved in the collapsing 
spherical region, and the virial radius R viT can be computed relatively simply by following 
the kinetic and potential energies of the collapsing matter and the potential energy of the 
dark energy contained in the spherical region. This case of dark energy not participating 
in the virialisation process was examined in e.g. [16, 19]. However, when clustering dark 
energy also partakes in the virialisation process, the condition for virialisation changes 
from these simple prescriptions, because unless c 2 s = 0, dark energy is not conserved in 
the spherical region. This necessitates a reexamination of the virialisation conditions 
in clustering dark energy cosmologies. Such an analysis has been carried out in [16], 
where both the homogeneous background as well as the clustered part of dark energy 
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take part in virialisation. See also [20-23] for other studies of dark energy virialisation. 
Based on our definition of Mq in equation (2.12), however, we will assume that only 
the clustered part of dark energy virialises. 

Our starting point is the tensor virial theorem for a system of point particles [24], 



ld 2 I 



jk 



2 (jf2 2K jk + W jk , (2.15) 

where Kjk and Wjk are the kinetic-energy and potential-energy tensor respectively, and 
Ijk is defined as 

Ijk = ^ ^ TTlg-Egj-Eaki (2'16) 
a 

with a running over the particles, and x a denoting their positions. Taking the trace 
of the equation (2.15) and assuming that the system is in a steady state, we find the 
well-known scalar virial theorem, 
d 2 I 

U-°< <2 - 17) 

where / = trace (I). Since the trace I is precisely the moment of inertia of the system, 
the condition for virialisation can also be paraphrased as the vanishing of its second 
time-derivative. 

For a spherical system, the moment of inertia for a set of non-relativistic particles 

is 

I = 2M tot R 2 /5, (2.18) 

where M tot denotes the mass that takes part in virialisation. Thus, in the case that mass 
is conserved in the spherical region (i.e., dM tot /dt = 0), equation (2.17) together with 
equation (2.18) provide a simple relation between the radius of the system and the first 
and second derivatives of the radius with respect to time at the time of virialisation: 

1 fdR\ 2 ld 2 R , 

In an Einstein-de Sitter universe (fl m = 1) where only dark matter clusters and virialises, 
mass is conserved in the collapsing region, and equation (2.19) dictates that virialisation 
happens when the radius of the system reaches half the radius at turn-around, i.e., 
i? v ; r = (l/2)i? ta , where turn-around refers to the instant at which the expansion of the 
spherical object stalls and the region starts to collapse (i.e., dR/dt\ tll = 0). 

Generalising to the case of mass non- conservation (such as when dark energy 
also clusters and virialises), we find from equations (2.17) and (2.18) the virialisation 
condition 

1 d 2 M tot | 2 dM tot dR | 1 UR\ 2 | ld 2 R = Q (22Q) 



2M tot dt 2 M tot R dt dt R 2 \ dt J R dt 2 
where M tot is the total mass of the system given by equation (2.14). In practice, after 
solving the spherical collapse, we locate the time at which equation (2.20) is satisfied. 
The total mass of the system at that time is then taken to be the virial mass M vir of 
the collapsed object. 
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Figure 1. The virial radius i? V ir, normalised to the radius at turnaround i?ta, as a 
function of the dark matter halo mass Mhaio- The left and right panels correspond to 
w = —0.8 and w = —0.9 respectively, while the upper and lower panels correspond 
to collapse at z = and at z = 1. The solid black lines correspond to (? s = 10 -6 , 
dashed lines to 10 -5 , dash-dot to 10 -4 , and dash-dot-dot-dot to 10 -2 . The solid red 
line corresponds to c~ = 1. 



2.5. Results 

Figure 1 shows the radius of the spherical halo at the time of virialisation as a function of 
the halo's dark matter mass for several clustering dark energy cosmologies. For clarity 
we have normalised the results to the radius at turn-around (i.e., the time at which 
dR/dt = 0). As with the linear threshold density S c in figure 5 of [9], a distinct step-like 
feature at the Jeans mass can be seen in the mass-dependence of the virial radius i? vir . 
The case of c 2 s = 10~ 6 , for example, has Mj ~ 10 14 M , as is evident in figure 1. 

The results presented in figure 1 can be understood in terms of a combination 
of several different effects: (i) the overall effect of the background homogeneous dark 
energy, (ii) the evolution of the ratio of the dark energy to dark matter densities, and 
(iii) the effect of the dark energy clustering. 
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Firstly, adding a homogeneous dark energy to the universe adds a negative pressure 
and hence a negative contribution to the spherical system's potential energy, making 
it more difficult for the system to virialise. The equation of state of dark energy tells 
us how strong this pressure is. Thus, the overall effect of the background dark energy 
is that Ryir/Rta becomes smaller compared to the Einstein-de Sitter case, in which 
i? vir /i?ta = 0.5 as discussed in the previous section. 

Secondly, the more w deviates from —1 (in the positive direction), the less rapidly 
the ratio Pq/ p m evolves with time. Fixing the ratio at z — 0, this means that at high 
redshifts dark energy is present in smaller amounts in models with w = —0.9 than in 
models with w = —0.8. The net effect is that R v i r /R t3u changes more dramatically as a 
function of time for w = —0.9 than for w = —0.8. (Note that the redshift at which dark 
energy starts to dominate over dark matter is z ~ 0.45 and z ~ 0.52 for w = —0.9 and 
w = —0.8 respectively.) 

The final effect of dark energy clustering draws i? vir /i? ta more closely to the value 
for an Einstein-de Sitter universe, because a larger portion of the spherical object now 
participate in virialisation, and the virialisation of dark energy has been assumed in our 
work to proceed like non-relativistic matter. As previously explained, the dark energy 
sound speed sets the scale at which dark energy clustering becomes relevant, and the 
transition from weak to strong clustering is clearly evident in figure 1 in the case of 
c 2 s = 10~ 6 . The accompanying figure 2 shows the contribution of dark energy to the total 
halo mass relative to the dark matter contribution at the time of virialisation. Clearly, 
dark energy clustering becomes more efficient with increasing halo mass, in the sense 
that the dark energy contribution to the total mass of the halo becomes increasingly 
important. However, this relative contribution does not increase indefinitely: Mq/M hll \ 
saturates at just above the Jeans mass, reaching a value of order 10" 3 . Comparing the 
Mg/M halo ratio at different redshifts, we see that the contribution of dark energy to the 
mass of the system is larger at z = than at z = 1. This follows simply from the fact 
that dark energy dominates the energy budget more at later times. 

3. The cluster mass function 

The cluster mass function is a power probe of dark energy properties. Firstly, its redshift 
dependence can be used to pin down the time evolution of the linear growth function, 
and hence the also the dark energy density as a function of time. Secondly, because 
the cluster mass function counts the comoving number density of virialised objects per 
mass interval, it can also be used to probe scale-dependent effects, such as the presence 
of a sound horizon due to a sound speed in the dark energy fluid. 

A fully accurate cluster mass function can in principle only be computed by way of 
numerical simulations (e.g., iV-body methods). However, for the purpose of a parameter 
error forecast, it suffices to estimate of the generic effects of dark energy using semi- 
analytical methods. To this end, we adopt the Press-Schechter formalism [14], which 
prescribes a cluster mass function in the mass interval [M vir , M viT + dM vir ] and at redshift 
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Figure 2. Dark energy contribution to the total mass at virialisation relative to the 
dark matter mass. The left and right panels correspond to w — —0.8 and w = —0.9 
respectively, while the upper and lower panels correspond to collapse at z = and 
z = 1. The solid black lines correspond to c 2 . — 10~ 6 , dashed lines to 10~ 5 , dash-dot 
to 10 -4 , and dash-dot-dot-dot to 10 -2 . The solid red line corresponds to c 2 . = 1. 



z of the form 



dn 



(M vir ,z)dM v 



2 Pm,C 



din v 



exp 



v 

T 



dM v 



(3-1) 



dM vir v "'" vll ''" / """ v " VvrM v 2 ir d\nM vlT 

Here, p mj0 is the present-day matter density, and v = 5 c {z) / a m {M Viri z), with 8 c (z) the 
linear threshold density (for matter) at collapse and a m (M v i r , z) the variance of the linear 
matter density field. The latter can be computed from linear matter power spectrum 
P^(k,z) via 

1 f°° 

a 2 m (M VlT ,z) = — dkk 2 \W(kR)\ 2 P l ^(k,z), (3.2) 

where |W(fcii)| is the Fourier space representation of a spherical top-hat filter of radius 
R, and M vir = ^-p m (z)R 3 relates R to the virial mass of the collapsed object via the 
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mean matter density p m (z) at redshift z.% The linear power spectrum P^(k,z) can 
be obtained from a Boltzmann code such as CLASS [15]; the publicly available version 
already incorporates an option for constant adiabatic dark energy sounds speeds. 

Apart from its effect on the linear power spectrum, a dark energy sound speed 
also gives rise to a mass-dependence in the linear threshold density 8 C [9]. This can be 
incorporated into the Press-Schechter mass function (3.1) simply by promoting 



where S c (M viT , z) is computed from following a spherical collapse as discussed in section 2. 
See [9] for details. Note that the definition of M vir here includes the contribution from 
clustered dark energy Mq. 

3.1. Results 

Figure 3 shows the cluster mass functions at redshifts z = and z = 1, for two 
cosmologies with w = —0.9 and w = —0.8 respectively, and a common sound speed 
c 2 — 1. Other cosmo logical parameters have been fixed at the WMAP7 ACDM best-fit 
values [1] . The corresponding linear matter power spectrum are also shown in the same 
figure. Dark energy dominates earlier in the case of w — —0.8 compared to w — —0.9, 
leading to a slightly scale-dependent suppression in the matter power spectrum of order 
8% at z = and 5% at z — 1. This suppression of power due to a larger w is also 
reflected in the cluster mass function in the form of a decrease in the number of massive 
objects. At 1O 16 M , the w = —0.9 cosmology has four times more objects than the 
w = —0.8 case at z — 0, and two times at z — 1. 

Compared with the effect of w, changing the dark energy sound speed from c 2 — 1 
to a smaller value has a considerably smaller effect on both the linear matter power 
spectrum and the cluster mass function. Figure 4 shows the linear matter power spectra 
for fixed w = —0.9 and w = —0.8, while varying the dark energy sound speed between 
c 2 = 10 -2 and 10~ 6 . The normalised (to the c 2 = 1 case) spectra show that changing 
the dark energy sound speed induces an additional scale dependence the matter power 
spectrum. In particular, on scales above the sound horizon and below the Hubble 
horizon, the amplitude of the matter power spectrum is enhanced. In the w = —0.8 
case, the enhancement is only 2% at z = 0, and less than 1% at z = 1. The w = —0.9 
case sees even less enhancement. Note that on scales close to the Hubble horizon (i.e., at 
very small values of k), a c 2 s <C 1 in fact causes the matter density contrast to decrease. 
This is an artifact of the synchronous gauge used in the computation of the matter 
power spectrum (see, e.g., [10]), and has no consequence on the observable scales. 

Interestingly, the effect of a dark energy sound speed appears to be somewhat more 
sizable on the cluster mass function. Figure 5 shows the cluster mass functions for 
a range of dark energy sound speeds, all of which have been normalised to the c 2 s = 1 

| On average, the dark energy does not contribute to the mass, since only the perturbations have been 
included in the definition Mq in equation (2.12). 




cr m (M vir ,z) 



a m (M vil , z) ' 



(3.3) 



Confronting the sound speed of dark energy with future cluster surveys 



11 




Figure 3. Cluster mass functions (left) and linear matter power spectra (right) for 
the case of c 2 s = 1. Solid lines represent a cosmology with w = —0.8, and dashed lines 
to w = —0.9. The upper and lower panels correspond to collapse at z = and z = 1 
respectively. 



case. Evidently, clustering dark energy with sound speeds in the range c 2 s = 1CT 6 — > 10~ 4 
increases the number of massive objects (~ 1O 16 M ) at z = by as much as 15% in the 
w = —0.8 case, and 5% in the w = —0.9 case. 

The astute reader will have noticed by now that some of the cluster mass functions 
presented in figure 5 have rather peculiar shapes. This calls for an explanation. Firstly, 
it is clear from equation (3.1) that the Press-Schechter mass function represents an 
interplay between the linear threshold density at collapse 5 C and the variance of the 
matter density field a m , with the relevant measure being the ratio v = S c /a m which 
generally increases with M vir . The shape of the cluster mass function for large masses is 
dominated by the exponential factor exp(— z/ 2 /2) in equation (3.1), as can be discerned 
from figure 3. Here, for a fixed value of M vir , the role of dark energy clustering is 
to subtly increases <7 m (M vir ) via its effect on the linear matter power spectrum as per 
equation (3.2), especially for cluster masses exceeding the Jeans mass Mj. This leads to 



Confronting the sound speed of dark energy with future cluster surveys 



12 




1 .00 



.ft °- 98 

Q_ 



0.96 



0.94 



10" 4 10" 



^o~ 2 10" 

k/Mpc" 1 



10 L 



10" 4 10" 



j i i 

2 i n-1 



10"^ 10 
k/Mpc"' 



10 L 



Figure 4. Linear matter power spectra for various values of (? s normalised to the = 1 
case. The left and right panels correspond to w — —0.8 and w = —0.9 respectively, 
while the upper panels show 2 = and the lower ones z = 1. The solid black lines 
correspond to c 2 = 10~ 6 , dashed lines to 10 -5 , dash-dot to 10 -4 , and dash-dot-dot-dot 
to 10~ 2 . 



a lower v, and hence an increased number of halos on large scales. This is the primary 
cause of the shapes of the cluster mass functions in figure 5. 

Secondly, when the dark energy sound speed deviates from c\ — 1, the linear 
threshold density S c is a function of the cluster mass M vir . This dependence on M vir is 
most pronounced at masses around the Jeans mass Mj\ away from Mj, 5 C asymptotes to 
5f n at M vir < Mj, and to 5™ ax at M vir > Mj [9]. For the case of c 2 s = 10" 6 (solid line 
in figure 5), for instance, the sudden increase in S c occurs at around Mj = 4 x 1O 14 M . 
This increase in 5 C compensates the growth of a m with respect to M vir , thereby causing 
a plateau in the cluster mass function. At M vir ^> Mj, 5 C asymptotes to £™ ax , and the 
cluster mass function recovers its exponential behaviour. 

To clarify this interplay between <5 c (M vir ) and cr m (M vir ), we show in figure 6 the 
cluster mass functions that would have arisen had the linear threshold density 5 C been 
fixed at 5™ (i.e., the value of 8 C at in the case of no dark energy clustering). Observe 
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Figure 5. Cluster mass functions for various values of <? s normalised to the <r s = 1 
case. The left and right panels correspond to w = —0.8 and w = —0.9 respectively, 
while the upper panels show z — and the lower ones z = 1. The solid black lines 
correspond to c 2 s — 10~ 6 , dashed lines to 10 -5 , dash-dot to 10 -4 , and dash-dot-dot-dot 
to 10~ 2 . 



the absence of the aforementioned plateau. Note also that the number of very massive 
clusters is a little higher here than in the figure 5. This is because the higher 5 C at large 
M vir in figure 5 leads to a stronger exponential suppression. 

Finally, we note that the mass-dependence of the linear threshold density 5 c (M vir ) 
has been computed in this work using a linear treatment of the dark energy perturbations 
adopted from [9]. While this approximation works well up to the Jeans scale, it breaks 
down at M > Mj because the dark energy density contrast inevitably exceeds unity at 
the time of collapse. This break-down leads to an overestimation of the exact value of 
S c at M > Mj, because nonlinear dark energy clustering is expected to feed back on the 
matter clustering more efficiently, thereby leading to an earlier collapse for w > — 1 than 
is implied by our linear treatment. An earlier collapse in turn leads to less evolution for 
the linear matter density contrast, and hence a 5 C at M > Mj that is lower than that 
computed using our linear approximation. 
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Figure 6. Same as figure 5, but for a linear threshold density 8 C fixed at a mass- 
independent S c (z) — 6 c (z, Cg = 1). 



Fortunately, however, the mass-dependence of 5 C works to counter the effect of dark 
energy clustering on a m , so that, as shown in figures 5 and 6, the contrasts between the 
cluster mass functions for different values of c 2 s are smaller when S c has been allowed to 
vary with M vir than when it has not. (Compare figures 5 and 6 also to figure 7 of [12].) 
In other words, our overestimation of 5 C at M > Mj leads to an underestimation of the 
effect of (? s on the cluster mass function. Thus, when we confront our semi-analytical 
cluster mass functions with future observations in the next section, the sensitivities to 
c\ that will be derived will be conservative. 

4. Constraining dark energy with future cluster surveys 

In this section we perform a forecast for the sensitivity of the Euclid cluster survey to 
the dark energy equation of state parameter w and non-adiabatic sound speed c 2 s using 
the Fisher matrix formalism. As we expect data of the cosmic microwave background 
anisotropies from Planck to be available long before Euclid is launched, we combine the 
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constraints from the Euclid cluster survey with those from Planck. 



4-1. Fisher matrix formalism 

The Fisher matrix formalism provides a quick way to compute the expected sensitivity of 
an experiment to cosmological parameters. Although Monte Carlo methods are arguably 
a more robust forecast tool [25] , we opt not to use it here because of the time consuming 
real-time evaluation of the spherical collapse. Furthermore, the Fisher matrix has the 
advantage that, once computed, individual Fisher matrices for different observational 
probes can be combined together by simple addition. Priors on specific parameters can 
likewise be incorporated into the Fisher matrix by adding cr~ 2 (pi) to the ii entry. We 
shall make use of these properties when assessing the constraining power of Euclid and 
Planck, both separately and in combination. 
The Fisher matrix is defined as 

where C is the likelihood function, and Pi is the ith parameter of the model. The 
inverse Fisher matrix, (F^ 1 )^, provides the best attainable covariance matrix, regardless 
of the specific method used to estimate the parameters from the data [26]. As a 
consequence, (F" 1 )^ 2 is the optimal statistical uncertainty attainable on parameter 
Pi after marginalisation over all other model parameters. 

Assuming a Gaussian likelihood function, the Fisher matrix for the Euclid cluster 
survey can be constructed as [27, 28] 

F .._y^^J_ (42) 



where 



dpi dpj N a ' 

d 2 V f°° dn(M,z, 



N n = AttAz 



dVLdz 



is the number of clusters with masses above the detection threshold M min (z a ) in a 
redshift bin centered on z a . In our analysis we consider 13 redshift bins from z — 0.1 to 
z = 1.4 with bin width Az = 0.1, chosen so that Az is twice the expected uncertainty 
in the determination of the cluster redshifts in the survey [17]. We assume the redshift 
bins to be top-hat functions, effectively ignoring small leakages that may arise as a 
result of a cluster being assigned to the wrong redshift bin. Note that a cluster survey 
such as Euclid extends to z ~ 2. However, our analysis shows that the number of 
detectable clusters is negligible already at z ~ 1.4 since M min at this redshift is of 
order 1O 16 M . The remaining quantities in the expression (4.3) are AQ the solid angle 
covered by the survey, d 2 V/ (dQdz) (z a ) the comoving volume element at redshift z a , 
and dn/dM(M,z a ) is the cluster mass function discussed in section 3. 
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Table 1. Euclid survey parameters. See section 4.2 for an explanation. 



Parameter 


Value 


Parameter 


Value 


An 


20000 


deg 2 


Cnfw 


5 


G 


1 


arcmin 


a 


2 


n bg 


30 


arcmin -2 




2 




0.1 






1 


Az 


0.1 
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4-2. Euclid survey parameters 

In order to determine the detection threshold of the weak lensing survey, i.e., M min (z a ) 
in equation (4.3), we adopt the approach of [27] and [29]. The average shear within a 
Gaussian filter of angular size 9 G) is related to the mass of the cluster M vir by [27], 

M vir / (irR 2 ) 



k g = a (9 G ) 



(4.4) 



Assuming the mass distribution in each cluster obeys a Navarro-Frenk- White density 
profile with a scale radius R s = i? vir (M vir )/c n f w , where i? vir (M vir ) is determined in this 
work via the spherical collapse, the coefficient a(9 G ) is given by 

J °° dx (x/x 2 G ) exp (-x 2 /x 2 G ) f (x) 



a(9 G ) 



(4.5) 



In (1 + c nfw ) - c nfw / (1 + c nfw ) 

Here, x = 9/9 s is an angular coordinate scaled with the angular scale radius 9 S = 
R s /dA (zi), dA (zi) the angular diameter distance to the cluster at redshift zi, x G = 9 G /9 S 
the smoothing scale, and the dimensionless surface density profile / (x) is given by 
equation (7) of [29]. 

In a geometrically flat universe, the mean inverse critical surface mass density is 



s- 1 = 



4ttG 



-X (zi) 



P° dz dn/dz (1 - x (zi) /x(z)) 



(4.6) 



C (1 + zi) - n bg 

where x denotes the comoving radial distance to the cluster, and dn/dz is the number 
density of source galaxies per steradian at redshift z, normalised such that = 
J °° dz dn/dz. We assume that dn/dz takes the form [30] 



dn 
dz 



dz = n h „ 



z T \ z ° 



exp 



dz, 



(4.7) 



where 7, (3, and z are chosen to fit the observed galaxy redshift distribution. In our 
analysis, we fix the parameters at 7 = 2, (3 = 2, and z = 1. 

In order for Euclid to detect a lensing signal, the shear k g must be larger than the 
noise. Here we adopt a detection threshold of k g = 4.5<7 no i se , where the noise term o" no i S e 
is modeled as a ratio of the mean ellipticity dispersion of galaxies a e and the number of 
background galaxies within the smoothing aperture N hg = 4irn hg 9 G [31], i.e., 
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Table 2. Fiducial parameter values of the two models analysed in this work. 



Parameter 


Model 1 Model 2 


w 


Dark energy equation of state 


-0.9 -0.8 


log c 2 s 


Dark energy sound speed 


-6 -6 


h 


Hubble parameter 


0.710 


ioon b h 2 


Physical baryon density 


2.26 




CDM density 


2.22 


log 10 10 A|, 


Primordial fluctuation amplitude 


1.39 


n s 


Scalar spectral index 


0.963 



All parameters relevant for the Euclid analysis are listed in table 1. Parameter values are 
sourced from primarily from [17]; parameter values missing in [17] are taken to match 
those for the LSST [27]. 

The detection threshold M min (z) at a given redshift z an now be computed from 
equation (4.4) simply by finding the lowest value of M vir at z that satisfies the condition 
k g > 4.5<T no i SC . For cosmologies in the vicinity of ACDM, the detection threshold is 
M min ~ 1O 13 M at z = 0.15, and M min ~ 9 x 1O 15 M at z = 1.15. 

4-3. Fiducial Cosmologies 

We consider a 7-parameter cosmological model defined by the free parameters 

9 = {w, log c 2 s , h, Vl b h 2 , n cdm , log 10 10 A^, n s }, (4.9) 

denoting, respectively, the dark energy equation of state, its sound speed, the Hubble 
parameter, the physical baryon density, the CDM density, the scalar fluctuation 
amplitude, and scalar spectral index. A flat spatial geometry, i.e., Qq = 1 — Q c dm — 
has been assumed. We examine two different fiducial cosmologies in our analysis: 

(i) model 1: w = -0.9 and c 2 s = 10~ 6 , 

(ii) model 2: w = -0.8, and c 2 s = 10~ 6 , 

where the remaining model parameters have fiducial values set by the ACDM best-fit 
of the WMAP 7-year observations [1]. See table 2 for details. The value of c 2 = 10 -6 
has been chosen such that the characteristic sound horizon scale falls squarely within 
the length scales resolved by the Euclid cluster survey (recall that the Jeans mass for 
c 2 = 10" 6 is ~ 1O 14 M ), while w = —0.9, —0.8 correspond approximately to the la and 
2a lower limits inferred from the WMAP7+BAO+ifo observations [1]. 

Ideally, we would also like to examine Euclid's sensitivity to dark energy models 
with sound speeds close to unity. However, we find that this exercise is not feasible 
within the Fisher matrix framework. This is because, as shown in figures 5 and 6, a 
sound speed c 2 larger than ~ 10~ 2 has essentially no observable effect on the cluster 
mass function on the mass scales probed by the Euclid cluster survey (1O 13 M < M < 
1O 16 M ). The corresponding likelihood function is virtually flat in the c 2 direction in 
the neighbourhood of the fiducial model, which is ultimately translated by the Fisher 
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Table 3. la sensitivities of the Euclid cluster survey to various cosmological 
parameters in the two fiducial models considered in this work. Note that the 
uncertainties in lOOOfc/i 2 and n s are identically the priors we set on these parameters 
using the WMAP 7-year data. 



Parameter 


Model 1 


Model 2 


w 


0.054 


0.065 


log c 2 s 


1.19 


0.46 


h 


0.074 


0.089 


ioon b h 2 


0.10 


0.10 




0.039 


0.038 


log 10 10 A|, 


0.12 


0.14 


n s 


0.020 


0.020 



matrix into unbounded or unrealistically lax constraints on c 2 s (recall that a Fisher 
matrix forecast is based upon evaluating the curvature of the likelihood function locally 
at the fiducial model). The failure of the formalism to take into account the full extent 
of the likelihood function is a severe drawback of the Fisher matrix approach, § which 
can be remedied only with a full Monte Carlo simulation [25]. We defer this exercise to 
a later study. 

4-4- Constraints from Euclid alone 

Table 3 shows the la sensitivities to various cosmological parameters attainable by the 
Euclid cluster survey alone. Note that when deriving these numbers, we have included 
the priors a (Vl b h 2 ) = 0.001 and a (n s ) = 0.02 on VL b h 2 and n s respectively, because the 
Euclid cluster survey has no particular constraining power for these parameters. These 
prior ranges correspond approximately to twice the la uncertainties inferred from the 
WMAP 7-year data [1]. 

The Euclid cluster survey alone is already fairly sensitive to the dark energy 
equation of state and hence the evolution of the dark energy density, with a la sensitivity 
to w of better than 10% in both models 1 and 2. These sensitivities should be compared 
with the uncertainties on w from the current generation of measurements: for a flat 
model with c 2 s — 1, a constant w parameter can be constrained to the same level 
of uncertainty (a(w) ~ 0.05) only with a combination of the WMAP 7-year data, 
measurements of the baryon acoustic oscillations, and type la supernovae; dropping the 
supernova measurements immediately enlarges the uncertainty by at least a factor of 
two [1]. 

Our results also show that, if the fiducial dark energy sound speed falls well within 
the observable range of the cluster survey, then the model can be distinguished from 
c 2 — 0,1 with a good degree of certainty. In our example of logc^ = —6, the la 

§ Although the Fisher matrix approach is local in nature, it does return reasonable estimates if the 
likelihood function does not deviate too much from a (multivariate) Gaussian distribution, for which 
the curvature in a chosen direction is constant no matter the distance from the best-fit. 
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Table 4. Survey specifications of Planck in three frequency channels [35]. 





fOO GHz 143 GHz 217 GHz 


9 C (arcmin) 
CT T ,c (a» k ) 
<7 E ,c (/"K) 

4 


10.7 8.0 5.5 
5.4 6.0 13.1 
11.4 26.7 
757 1012 1472 



sensitivity of Euclid to \ogc 2 s is 1.19 and 0.46 for fiducial w values of —0.9 and —0.8 
respectively: the sensitivity to log c 2 s is better in the latter model because dark energy 
clustering is generally suppressed by a factor (1 + w) relative to dark matter clustering. 



4-5. Constraints from Planck+Euclid 



Planck measures the temperature and polarisation anisotropies of the Cosmic Microwave 
Background (CMB). These measurements are quantified at the simplest level in terms 
of auto- and cross-correlations of the temperature T and E-fB-tjpe polarisation 
fluctuations. Neglecting B-mode polarisation, the CMB Fisher matrix is given by [32] 



EE 

£ X,Y 



dpi 



dpj ' 



(4.10) 



where X and Y runs over TT, EE and TE, and the entries in the symmetric covariance 
matrix are given by 



Cov (CtT,£, Ctt,£ 
Cov (Cee,£, Cee,, 
Cov (Cte,£, Cte, 

Cov (Ctt,^ Cee,£ 
Cov (Ctt,£, Cte,£ 
Cov (Cee,£, Cte,£ 



(C T t,£ + B T 2 e ) , 



(2£+l)/ sky 
2 

(2^+1) /sky 

1 

(2£+l)/sky 

2 



{Cee,£ + B E 2 £ ) 2 , 

[ C TE,£ + (CTT,£ + B^ 2 ) (CeE,£ + B^ 2 )] , 



(21 + I) U y C ^ 



(2£+l)f sky 
2 



Cte,£ (Ctt,£ + B T fy , 
Cte,£ (Cee,^ + B E 2 ^ . 



(4.11) 



(2£+l)f sky 

Here, / s k y is the fraction of the sky remaining after removal of foregrounds and galactic 
plane contaminations, etc., and B x ,e denotes the expected error of the experimental 
apparatus [33, 34], 

4,^E(^)" 2 e^ +1)/ ^ (4-12) 

c 

where the index c labels the different frequency channels, o~ Xc is the variance of the 
instrumental noise in the temperature/polarisation measurement, 9 C is the width of the 
beam assuming a Gaussian profile, and i c = 2v^2 ln2/^ c is the corresponding cut-off 
multipole. Table 4 lists the values of these quantities specific to Planck's observations. 
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Table 5. la sensitivities to various cosmological parameters from the Planck 
temperature and polarisation measurements (TT, EE and TE). 



Parameter 


Model 1 Model 2 


a (w) 

cr(l0g C 2) 

a(h) 

a (ioon b h 2 ) 
<t (ion cdm ) 

a(lo g 10 10 A^) 

a(n s ) 


0.21 0.21 
5340.77 948.04 
0.082 0.090 
0.026 0.025 

0.52 0.57 
0.010 0.010 
0.0067 0.0063 



Table 6. la sensitivities to various cosmological parameters from combining the Euclid 
cluster survey and the Planck temperature and polarisation measurements (TT, EE 
and TE). 



Parameter 


Model 1 


Model 2 


w 


0.0060 


0.0046 


log c 2 s 


0.58 


0.22 


h 


0.0025 


0.0022 


ioon b h 2 


0.017 


0.017 




0.016 


0.014 


log 10 10 A|, 


0.0047 


0.0047 


n s 


0.0032 


0.0032 



Table 5 shows the la sensitivities of Planck to the parameters of our two fiducial 
models. Compared with the Euclid-only results in table 3, it is immediately clear 
that Planck has very little sensitivity to dark energy parameters compared with 
cluster surveys; the la sensitivity to w is only about 20%, while logc^ is completely 
unconstrained, consistent with previous expectations (e.g., [36-38]). Planck's sensitivity 
to h is comparable to that of the Euclid cluster survey (~ 10%), but the latter performs 
ten times better for the fl c dm measurement. 

What is more interesting is when observations from Planck and the Euclid cluster 
survey are used in combination. The combined la sensitivities of Planck and the Euclid 
cluster survey are displayed in table 6. Compared with the expectations from Planck or 
Euclid alone, we see that Planck+Euclid improves the sensitivities to some parameters, 
especially h and w, by more than tenfold; the la sensitivities to w and h from the 
combined analysis are both now better than 1%. Such a significant improvement 
shows that Planck and the Euclid cluster survey are individually sensitive to different 
degenerate combinations of these parameters. However, once these observations are 
used in combination, the degeneracies are completely broken. Figure 7 illustrates how 
the degeneracies in h and w inherent in the Planck and in the Euclid observations can 
be broken by a combined analysis. 

Importantly, while Planck alone has no constraining power whatsoever on the dark 
energy sound speed, because of the breaking of parameter degeneracies discussed above, 
the combination of Planck and Euclid turns out to improve the Euclid-alone sensitivity 
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0.010 



0.005 



§ 0.000 



-0.005 



-0.010 

-0.04 -0.02 0.00 0.02 0.04 
Aw 

Figure 7. lcr sensitivities in the (h, w)-parameter space for the Euclid cluster survey 
(black lines), Planck (red), and combined Planck+Euclid (blue), assuming fiducial 
model 2. The dotted lines represent sensitivities after marginalising over all other 
cosmological parameters, while the solid lines correspond to no marginalisation (i.e., 
parameters not shown in this plot are kept fixed at their fiducial values). 

to log c 2 by about a factor of two. Parameters such as the physical baryon density VL^h 2 
and the scalar spectral index n s , to which Planck is most sensitive, also benefit somewhat 
from a combined analysis with Euclid: the sensitivity to n s from the combined improves 
on the Planck-alone result by about a factor two, while for Q^h 2 we find a ~ 30% 
improvement. 

5. Discussions and conclusions 

Current observation of an accelerated universal expansion is commonly attributed to 
the presence of a dark energy component in universe's total energy budget. However, 
apart from that it should have some negative pressure, our understanding of the details 
of this dark energy component is quite limited. In this work, we have investigated how 
future cluster surveys can help to shed light on this dark energy through their potential 
to observe galaxy clusters numbering in the hundred thousands. We have considered in 
particular how the cluster survey of the ESA Euclid project will allow us to determine 
the equation of state and the sound speed of a generic dark energy fluid. 

The main role of a dark energy sound speed is that it enables the clustering 
of dark energy on length scales above the sound horizon, thereby allowing the dark 
energy component to participate actively in the formation of structure in the universe. 
This means that, if the sound horizon is comparable to or smaller than the typical 
length scales of a galaxy or galaxy cluster, one should expect to find dark energy 
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bound gravitationally in the collapsed objects and contribute to their mass. We have 
demonstrated in this work, using the spherical collapse model, that the maximum dark 
energy contribution to the total cluster mass is of order a few tenths of a per cent at 
the time of virialisation. 

Secondly, because dark energy clustering feeds back on the evolution of the dark 
matter density perturbations, it leads to a scale-dependence in the dark matter clustering 
that turns out to be more significant than is suggested by the mere ~ 0.1% infall 
discussed immediately above. For dark energy sound speeds approximately in the range 
= 1CT 6 — > 1CT 4 , the effect of dark energy clustering is directly visible in the cluster 
mass function in the mass range 10 12 — > 1O 16 M ; for dark energy models consistent with 
present observations, the number massive clusters (~ 1O 16 M ) can potentially change 
by as much as 15% when compared with the no clustering case. 

Using our predictions for the cluster mass functions, we have performed a Fisher 
matrix forecast in order to assess the potential of the Euclid cluster survey to constrain 
dark energy parameters. We find that a la sensitivity to the (time-independent) 
dark energy energy equation of state parameter w at the sub-percent level is possible, 
especially when the the Euclid cluster survey is complemented by CMB observations 
from Planck. Furthermore, if the dark energy has a sound speed that deviates 
appreciably from unity, we will be able to pin it down with future cluster surveys to 
roughly within an order of magnitude. 

The results presented in this work are well in line with those from previous studies, 
where observational probes other than cluster surveys had been considered [7, 8, 10, 39, 
40]. We therefore conclude that future cluster surveys will be competitive with other 
probes in their constraining power on dark energy parameters, and that the optimal 
strategy is a combination of different observations. 
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